Shape of primary proton spectrum in multi-TeV region 
from data on vertical muon flux 
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It is shown, that primary proton spectrum, reconstructed from sea-level and underground data on 
muon spectrum with the use of QGSJET 01, QGSJET II, NEXUS 3.97 and SIBYLL 2.1 interaction 
models, demonstrates not only model-dependent intensity, but also model-dependent form. For 
correct reproduction of muon spectrum shape primary proton flux should have non-constant power 
index for all considered models, except SIBYLL 2.1, with break at energies around 10-15 TeV and 
value of exponent before break close to that obtained in ATIC-2 experiment. To validate presence 
of this break understanding of inclusive spectra behavior in fragmentation region in p-air collisions 
should be improved, but we show, that it is impossible to do on the basis of the existing experimental 
data on primary nuclei, atmospheric muon and hadron fluxes. 

PACS numbers: 95.85.Ry, 96.50.sb 



I. INTRODUCTION 

At present information on the characteristics of 
hadronic interactions in fragmentation region is still 
scarce or missing and experiments with 'roman pots' are 
anticipated to improve the situation. Some of this in- 
formation, in principle, could be obtained with the use 
of the data on cosmic ray (CR) muon and hadron spec- 
tra, provided primary spectra are known with high pre- 
cision, but that is not the case. The obvious obstacle 
here is that at high energies primary cosmic ray (PCR) 
fluxes, measured in direct experiments, themselves are 
functionals of various interaction parameters plus their 
accuracy is appreciably affected by additional systematic 
effects [il, 0, B 01 • In the series of our papers [3, [a] we 
underlined, that these effects can lead to underestima- 
tion of light nuclei fluxes and thus can explain discrep- 
ancy between measured and calculated muon fluxes for 
Ef, > 100 GeV. Prehminary data of ATIC-2 0, cover- 
ing the gap between magnetic spectrometer and emul- 
sion chamber experiments, seems to be in concordance 
with our conclusions, but situation is more complicated 
in fact, as further consideration will show. ATIC-2 slope 
of proton spectrum 7p = 2.63 [8] for primary energies 
below 10 TeV is in remarkable contradiction with pre- 
viously measured values 7^ = 2.74 and 7p = 2.80 by 
RUNJOB !9] and JACEE [l^l experiments correspond- 
ingly, but this discrepancy is removed by steepening in 
ATIC-2 proton spectrum at energies above 10 TeV. These 
new data were already exploited in extensive calculations 
of muon and hadron fluxes with number of interaction 
models at different atmospheric depths and zenith an- 
gles in [ll| where it was shown, that their use allows 



to get reasonable agreement with the most of the data 
under appropriate choice of hadronic interaction param- 
eters. In fact, it is not possible to obtain concordant 
conclusions on primary spectra or hadronic interactions 
parameters coming even from much smaller subset of the 
experimental data. In this paper we demonstrate this on 
the basis of analysis of the data on muon flux at vertical 
direction and only one set of the data on hadron flux of 
EAS-TOP [l2]. First, we use the data of sea-level and 
underground experiments to obtain conclusions on be- 
havior of muon spectrum at sea-level in the energy range 
40 GeV-10 TeV. Further we analyze influence of uncer- 
tainties in the muon data and interaction parameters on 
properties of reconstructed primary proton fluxes. And 
finally we show, that such different in approach and char- 
acteristics interaction models as SIBYLL 2.1 [13, llJ] and 
NEXUS 3.97 [l5| can bring to hardly distinguishable pre- 
dictions on muon and hadron fluxes. 



II. SEA-LEVEL MUON SPECTRUM FROM 
UNDERGROUND EXPERIMENTS DATA 

Depth-intensity relation, needed for reconstruction of 
sea-level muon spectrum, may be obtained via numerical 
solution of one-dimensional transport equation. In ad- 
joint approach this equation has the following form [iq 
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Here q{t, E) — is survival probability of muon with en- 
ergy E, being born at the distance t from detector, 
a — total interaction cross-section, Wp^E^E'), f3 — 
i, r,p,h — differential cross-sections for processes of ion- 



ization, bremsstrahlung, pair production and photonu- 
clear interaction correspondingly, D{t, E) — detector 
sensitivity function. The numerical method, applied for 
solution of this equation [17], allows to avoid any approx- 
imations (such as continuous losses one) and to obtain 
muon intensities at large depths of matter with account of 
fluctuations in all muon interaction processes. Accuracy 
of our calculations of muon survival probabilities and in- 
tensities was thoroughly examined in [l7[ and compari- 
son with the results of Monte-Carlo codes MUM '1^ and 
MUSIC [l^l is presented in Fig. [1] Anticipating further 
discussion of sea-level muon spectrum behavior it is nec- 
essary to note, that our calculations give upper estimate 
of muon flux at large depths in comparison with MUM, 
because of use of ~1% lower muon energy losses [l7|. 

To describe the data on muon intensity underground 
and at sea-level, and to estimate influence of uncertainty 
in muon flux data on reconstruction of primary proton 
flux, we used two parameterizations in the simple form, 
proposed in work of Bogdanova et al. [23|. Original fit 
for the vertical from [20] 

S^{E) = l%/{E+U^)/{E + 2.lf-\ (cm2.s-sr-GeV)-\ 

(1) 
provides good agreement with the data at sea-level 

(Fig. [2]), but leads to underestimation of the muon flux 
for the depths below 6 km w.e. (Fig. [3]). 

To match better underground data for the depths 2- 
6 km w.e. we shall also apply modified fit with slightly 
< 10% increased intensity in multi-TeV region 

S^{E) = 20.8/(£;+194.3)/^2-'^\ (cm^-s-sr-GeV)-^ (2) 

As it is seen from Fig. [31 for depths from 4 km w.e. up to 
8 km w.e., corresponding to ~ 2.5 — 10 TeV median muon 
energies at sea level, use of this spectrum provides good 
agreement with the data of LVD "2?!, BNO [22| and Fre- 
jus [23, [2J collaborations and leads to underestimation 
of data of MACRO [H and Soudan ^^ experiments. 
Our consideration will touch muon energies only above 
40 GeV, to exclude different effects, not related to high- 
energy hadronic interactions features, such as geomag- 
netic effect, influence of uncertainties in low-energy inter- 
action models, or even absorption of low energy muons 
in ground as in case of L3-I-C detector. 



III. CALCULATIONS TECHNIQUE OF 
ATMOSPHERIC MUON AND HADRON FLUXES 

Average numbers of hadrons A^h(-£'N,> -Eth) and 
muons Nf^{E-^,> Eth) with energies above Eth in EAS 
from primary nucleon of energy E-^ were obtained with 
the help of one-dimensional hybrid code CONEX [40, 
l41j in regime of cascade equations solution for in- 
teraction models QGSJET 01 [H, SIBYLL 2.1 [H, 
NEXUS 3.97 dl and QGSJET-II-03 [H,!!!,!!!. To get 
differential spectrum of hadrons (or in the same way of 
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FIG. 1 (color online). Muon intensity in water and in 
standard rock. Water (muon spectrum at sea level 
from [iMl): solid line — present work; closed 
circles — MUM |18|, error bars show uncertainty due to 
±1% variation of total muon energy losses. Standard 
rock (muon spectrum at sea level from [2^): dashed 
line — present work, open circles — MUSIC [29J . 



muons) for some energy Eq the following simple formula 
had been used: 

^h(^o) - [/h(> Eo - A£;) - /h(> Eo + AE)] /2AE. 



Here /h(> Eq) is the integral spectrum of hadrons 



(3) 



h{>Eo)= / SNiEn)Ni,iEN,> Eo)dEN 



for primary nucleon spectrum S'n(-Bn)- 

The interval width AE must be chosen to provide 
the difference between integral intensities in ([3]) to be 
much larger then calculations error. Test computations 
for AE = 0.01£;o,0.02£;o,0.05£;o and observation levels 
820 and 1030 g/cm^ for energy Eq = 1 TeV brought to 
differential flux values lying within 3% from each other. 
Change of upper integration limit Eoo in formula for 
/h(> -Eo) from lO" x Eq to 10^ x Eq gives less then 1% 
increase of differential flux. And, the last, increase of 
number of primary energy bins in 7Vh(i?N, > £-0) from 10 
to 20 per order introduces ^ 2% variation to differential 
flux value. The listed error sources partially compensate 
each other and the total error of our calculations does 
not exceed 3%. The calculation were performed for the 
set of energies, coinciding with the set from EAS-TOP 
experiment paper [l2l |. 



IV. SEA-LEVEL MUON FLUXES 

As a basic model of PCR nuclei spectra the parameter- 
izations from [36i] were chosen. Nuclei with A > A were 
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FIG. 2 (color online). Sea-level muon spectrum. Experimental data: |2|] MACRO 1995, [30] L3+C 2004, 
PI BESS-TeV, ^ Rastin, [2l| LVD 1998, [33^ Baksan 1992, [s^ ASD 1985, H MSU 1994, [2! Frejus 1994. Left: 
Solid line — muon spectrum ^, dashed line — muon spectrum ^. Right: Muon spectra at sea level for PCR 
spectra from [36[ with high hclimn flux. Muon spectrum parametrization ([T]) is shown with crosses. 
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FIG. 3 (color online). Vertical muon intensity in 
standard rock. Experimental data: [22| Baksan. 
[m LVD, [131 Crouch^ MACRoTlli, H^ Soudan 
1,2, [H] Boulby, fH, HjFrejus. Present work 
calculations: solid line — for muon spectrum ([T]), 
dashed line — for muon spectrum ([2]). Neutrino 
induced muon contribution is taken from !39l|. 



treated in the framework of the superposition model, high 
accuracy of this approach is well known and was checked 
by our calculations with the use of CONEX both for 
muons and hadrons once again. 

Comparison of the calculated muon fluxes with the ex- 
perimental data, presented in Fig. [21 reveal familiar pic- 
ture of high energy muon deficit. The reasons of its ap- 
pearance were considered in our previous papers [5|, [6[ 
and they still hold true regardless of the fact, that three 
more interaction models were included in our analysis. 
All interaction codes, except QGSJET 01, satisfactory 



describe data on muon flux only up to E^ '^100 GeV 
and then one by one fail to do it. Accounting that 
such muon energies correspond to primary energies above 
1 TeV, studied with balloon(satellite)-borne emulsion 
chambers, we related muon deficit to underestimation of 
primary light nuclei fluxes, taking place in these exper- 
iments [5[, [6| . Unfortunately, disagreement between the 
models in the muon fluxes also appears at energies around 
100 GeV, thus making impossible precise reconstruction 
of primary nucleon spectrum for Ep^m > 1 TeV. In fact, 
in such conditions there are no reasons to rule out any 
of the models, except QGSJET 01, which leads to re- 
markable disagreement with the experiment even in the 
range of reliable magnetic spectrometers data on PCR 
and muon spectra. 

To find out, why the models differ in the predicted 
muon fluxes let us consider quite characteristic energy 
1.29 TeV, where discrepancies between the models reach 
appreciable values and the data on muons from under- 
ground installations are yet quite reliable. Contributions 
of primary protons to the differential flux of muons of 
the given energy, presented in Fig. [H show, that spread 
in muon fluxes between the interaction models is entirely 
due to uncertainties in the description of tt , i^^-spectra 
in fragmentation region x = i?,r,A'/-E'prim > 0.1. Since in- 
clusive muon flux is sensitive nearly only to the character- 
istics of the very flrst primary particle interaction, hence, 
the harder these spectra are in the particular model, the 
larger muon intensity its use leads to. For the lower 
values of x, i.e. for ii^prim > 10 TeV, all the models 
give practically the same muon yields. As noted above, 
in view of uncertain situation with primary spectra for 
Eprim > 1 TeV, one can not give preference to any of 
the models in comparison with the others. If to demand 
the minimal disagreement with the direct measurements 
data on PCR spectra, then obviously SIBYLL 2.1 satis- 
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FIG. 4 (color online). Left: contribution of primary protons with energies i^prim to the niuon differential spectrum at 
sea level for _E^ — 1.29 TeV. Right: inclusive spectra p + A ^ tt^ + X and p + A ^ K^ + X (scaled down by 10) for 
incident proton with energy 10 TeV. 



ties this requirement the best, or, on the other hand, one 
could say that it provides the most acceptable descrip- 
tion of TT , K production spectra in p-air collisions in 
fragmentation region. We shall discuss this affirmation 
in detail below. 



V. PROTON SPECTRUM FROM MUON DATA 

From the previous consideration it is clear, that recon- 
structed fluxes of protons shall be higher than measured 
in emulsion chamber experiments, but can be compa- 
rable with recently obtained data of ATIC-2 group Q. 
We performed the reconstruction simply by picking up 
of appropriate primary proton flux parameters in order 
to minimize deviation of the obtained muon fluxes from 
the parameterizations H]) and ([2]). 

First, let us consider an attempt to minimize maximal 
deviation from spectra H]) , (HD for ^^ = 40 GeV- 10 TeV 
in assumption, that primary proton spectrum can be de- 
scribed by single power law function Jp — AE~'^p in the 
entire energy range 100 GeV-500 TeV. In Fig. [5] and in 
TableUone can find the results of this reconstruction, up- 
per and lower lines for each muon spectrum parametriza- 
tions in the Fig. [5] correspond to low and high helium flux 
fits [3y|. As expected, the obtained spectra are flatter, 
than measured by RUNJOB and JACEE groups, but, 
except for QGSJET 01 model, agree well with ATIC-2 
results Q. Figure [5] shows, that in this case it is possi- 
ble to achieve agreement with the fitted muon spectrum 
within 10%, but only SIBYLL 2.1 reproduces its shape 
correctly, with other models it is not possible to get right 
muon spectrum slope variation. We should also note, 
that this behavior turned out to be insensitive to the 
choice of helium flux parametrization. As a result, there 
is a dip in the ratio of the muon flux, obtained from fit- 
ting of proton spectrum, to the parameterizations H]), ^ 



for energies around 1 TeV, where the underground data 
were already underestimated (overestimation of muon 
flux at higher energies does not compensate this effect 
completely), and growth in small energies range, which 
brings to contradiction with low energy Ef^ < 40 GeV 
data. Parametrization Q was also projected to make 
muon spectrum slope variation less sharp, but this lead 
to problems with its fltting in ^ 100 GeV energy range, 
as it is seen from Fig. [5] (right panel) for SIBYLL 2.1 and 
NEXUS 3.97 models. The problem with muon spectrum 
shape matching is better illustrated by upper left panel 
of Fig. [71 where it is shown, that correct reproduction 
of muon spectrum for energies below 1 TeV with single 
power law proton spectrum leads to appreciable overesti- 
mation of muon flux at higher energies. There are three 
possible explanations or solutions of this problem. First, 
the discrepancy can be completely removed by choice of 
appropriate interaction parameters, e.g. similar to those 
in SIBYLL 2.1. Another argument, which can be given 
is that the data on muon flux for energies above 1 TeV 
are not so definite to claim their inconsistency with the 
calculations, but it does not look well supported by un- 
derground data (see Fig. [3] and calculations [ij, ISlf). 
And the last possibility is to assume, that primary pro- 
ton spectrum is not monotonous and either has sharp 
break or slowly changing exponent jp. Let us consider 
the latter assumption, which finds experimental [T] and 
theoretical @, [52, [53] justifications, in more detail. The 
results for the simple case with break (Fig. [7]), which 
allows to achieve correct description of muon spectrum 
shape with right asymptotic and deviation in flux value 
< 3%, show, that small difference between spectra ([T]) 
and ([2|) results not only in different proton intensities, but 
also in break positions. The latter lies for parametriza- 
tion ^ in the primary energy range 10-15 TeV, the 
change in power index reaches appreciable values up to 
Ajp w 0.15 for QGSJET 01 and QGSJET II models 
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FIG. 5 (color online). Primary proton spectrum. Experimental data: 46] Ryan et al., 31] BESS-TeV, [47] SOKOL, 
[a] RUNJOB, [M] MUBEE, [10.] JACEE, J9,] ATIC-2, [50] KASCADE (hadrons). Solid lines and dashed lines show 
primary proton spectra for muon flux parametrizations ^ and ^ correspondingly. Upper and lower lines reflect 
uncertainty in helium flux according to [36[ . 



(Table m. Proton spectra, obtained from muon flux ^, 
with QGSJET II and NEXUS 3.97 models are in the 
best agreement with ATIC-2 data, while SIBYLL 2.1 pro- 
vides intermediate between ATIC-2 and emulsion cham- 
bers experiments slope value. Spectra, reconstructed 
from paramctrization ([1]), have breaks at 3-6 TeV and 
in case of QGSJET II proton flux poorly agrees with 
experiments at primary energies around 100 GeV. Evi- 
dently, the latter problems are explained by too low, in 
comparison with underground data, muon flux and this 
parametrization is considered here mostly for estimation 
of sensitivity of primary spectrum features to the choice 
of reference muon flux. 

It is necessary to note, that due to low sensitivity of 
differential muon flux to helium and heavier groups of 
primary nuclei it is impossible to derive any conclusions 
on presence of the break in these PCR components. For 
illustration let us consider example of calculations for 
QGSJET 01 and high helium flux, where the break in 
proton spectrum is positioned at ii^br = 15 TeV and 
change of power index is equal to 0.14 (see Table |T|. In- 



troduction of rigidity-dependent break in He spectrum at 
2_Ebr per nucleus of the same value A7 = 0.14 gives re- 
markable discrepancy between calculated muon spectrum 
and parameterization ^ only for energies above 7 TeV, 
which reaches 10% at 20 TeV. To correct this asymp- 
totic behavior it suffices to reduce A7 to 0.11 simulta- 
neously for protons and helium without change of the 
break position, and thus we get proton spectrum lying 
well within corridor between parameterizations for high 
and low helium fits, shown in Fig. [7] Hence, this cor- 
ridor covers all possible cases of He flux behavior (with 
or without break), provided the helium flux stays within 
limits, given in ^3^]. 



Summarizing we can say, that primary proton spec- 
trum shape turns out to be sensitive to the choice of 
interaction model and allows presence of break at 10- 
15 TeV with Ajp up to 0.15, which can be shghtly soft- 
ened though, if to allow presence of the same break in 
other PCR components spectra. 
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FIG. 6 (color online). Ratios of muon fluxes, obtained 
for power-law primary proton spectra, given in Fig. [51 
to muon flux parameterizations ([l]) and ([2]). 



VI. HADRONS 

Comparison of our calculations of hadron spectrum for 
the primary spectra from [36| (high helium flux) with the 
recent measurements, performed by EAS-TOP collabo- 
ration [I4I, is presented in Fig. [5] (left panel). First, let 
us note the following facts. Below 100 GeV all calcu- 
lated spectra have breaks, caused by non-perfect match- 
ing of low-energy interaction model GHEISHA to the 
high-energy models. Shape of the measured hadron spec- 
tra also breaks at energies above 4 TeV and the data be- 
come less definite, thus in the forthcoming analysis we are 
going to consider the data for energies from 129 GeV to 
4 TeV. For these energies QGSJET 01, QGSJET II and 
SIBYLL 2.1 quite reasonably reproduce the shape of the 
measured hadron spectrum, NEXUS 3.97 leads to spec- 
trum with almost constant power index. One can see, 
that the most consistent description of the data for spec- 
ified energies provide QGSJET 01 and SIBYLL 2.1. In 
contrast with the muons there are no energy range, where 
the models agree on the hadron fluxes and the reasons 
of this disagreement are not as simply to point out as in 



TABLE I. Parameters of primary proton spectrum fits 
Jp — AE >", (cm^ • s • sr • GeV)^^, for muon spectrum 
paramctrization ([2]). 

Low fielium flux 



No break 



With break 



Model 



Below break Break Above break 

7 yl 7 (TeV) A 7 



QGSJET 01 
QGSJET II 
NEXUS 3.97 
SIBYLL 2.1 



1.42 2.675 0.93 2.620 

1.04 2.650 0.69 2.595 

1.22 2.680 1.01 2.655 

1.29 2.695 



15 3.25 2.750 
10 2.39 2.730 
13 1.78 2.715 



High helium flux 



No break 



With break 



Model 



Below break Break Above break 

7 yl 7 (TeV) A 7 



QGSJET 01 
QGSJET II 
NEXUS 3.97 
SIBYLL 2.1 



1.45 2.690 0.95 2.635 

1.08 2.670 0.69 2.610 

1.25 2.700 1.01 2.670 

1.37 2.720 



15 3.65 2.775 
14 3.18 2.770 
12 2.14 2.750 



the case with muons. The most important characteris- 
tics in this analysis are total inelastic cross section, de- 
termining chances of primary particle to survive, shapes 
of inclusive spectra p + A -^ p + X, p + A -^ n + X, 
TT* + A —^ IT + V in the very forward region, responsi- 
ble for substantial process of leading particles production 
(see Fig. [5] for the listed spectra) . Let us briefly outline 
the major conclusions, which one may come to in the 
given situation. NEXUS 3.97 gives the lowest fluxes as 
of hadrons in total, so of nucleons and mesons (Fig. [T0|) . 
and this happens in spite of the lowest inelastic cross- 
section values. Inclusive spectrum p + A —^p + X im- 
mediately helps to flgure out, that incident protons in 
NEXUS 3.97 have comparably low chances to save most 
of their energy in collision and this leads to such low nu- 
cleon flux, the same can be said about meson flux and 
production of pions by pions. Similarly, from comparison 
of the inclusive spectra, it can be easily understood, why 
QGSJET II gives the highest hadron flux. Note, that 
SIBYLL 2.1 concedes to QGSJET II in hadron inten- 
sity mostly because of less effective production of leading 
neutrons in p-a.ii collisions and due to the larger total 
interaction cross-section. 

Thus, from analysis of the data on hadron flux it is 
difficult to imply any strict constraints on inclusive spec- 
tra shapes, since mechanism of hadron spectrum forma- 
tion is more sophisticated than that in the case of muon 
spectrum. SIBYLL 2.1 and QGSJET 01 display quite 
different behaviors of the relevant inclusive spectra and 
total interaction cross-sections, but both models almost 
equally succeed in description of the EAS-TOP data (i.e. 
produce close hadron fiuxes). 
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FIG. 7 (color online). Upper left panel shows ratios of muon fluxes, obtained for primary proton spectra with 
constant power indices equal to those before break, to muon flux parametrization ([2]). Other designations as in 
Fig. El 
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FIG. 8 (color online). Hadron spectra at the EAS-TOP depth t = 820 g/cm^. Le/f; calculations for primary spectra 
from j36|| with high helium flux. Right: calculations for proton spectra from Fig. [7] (for SIBYLL 2.1 from Fig. [5]). 
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FIG. 9 (color online). Left: inclusive spectra n^+A^TT^+X,p + A^p + X (scaled down by 10), 

p + A ^>- n + X (scaled down by 50) for incident particles with energy 10 TeV. Right: contribution of primary 

particles with energies £^prim to the hadron differential spectrum Sh[Eh) at 820 g/cm^ for Eh — 1.29 TeV. 
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FIG. 10 (color online). Spectra of nucleons (left) and mesons (right) at 820 g/cm^ for primary spectra [36] with high 
helium flux. 



The given standard approach to analysis of situation, 
in fact, is of little sense, since it is based on assump- 
tion about validity of primary spectra in form of fits 
from [3y| , which was called into question in our previous 
discussion. In this case it is logical to analyze interaction 
models self-consistency, i.e. their ability to give correct 
estimates of several CR components at once. Provided 
we know behavior of primary proton spectra for every 
model, required to match the data on muon flux, we may 
check how these proton spectra agree with the data on 
hadrons. In Fig. [8] (right panel) we give hadron intensi- 
ties, calculated for primary proton spectra with breaks 
from Fig. [7] (for SIBYLL 2.1 see Fig. [51), corresponding 
to muon spectrum parametrization ^ . After increase of 
primary nucleon flux, dictated by the data on muons, 
one can see, that the best agreement with EAS-TOP 
measurement provide NEXUS 3.97 and SIBYLL 2.1. It 
would be interesting to note that two models with differ- 



ent philosophies and inclusive spectra give the most self- 
consistent results on muons and hadrons, but, of course, 
this conclusion must be taken with much care, since it 
is based on the single set of data and we have only indi- 
rect indications on the accuracy of this set, e.g. such as 
agreement of primary proton fluxes, obtained by EAS- 
TOP and KASCADE teams (the latter is derived from 
flux of unaccompanied hadrons [50|). If to try to perform 
the same analysis with the variety of the data, obtained 
at different atmospheric altitudes and zenith angles, no 
consistent notions of such kind will be obtained as can 
be easily seen from calculations, presented in 



11|. 



VII. CONCLUSIONS 



The progress in CR and high-energy physics, achieved 
during last 10-15 years allowed to turn from statements 



about satisfactory (qualitative) concordance between dif- 
ferent kinds of data to investigation of more fine ef- 
fects. As an example, we managed to show that recon- 
structed from the data on vertical muon flux primary 
proton spectra have not only expected interaction model 
dependent intensities, but also model-dependent shapes. 
It is demonstrated, that application of QGSJET 01, 
QGSJET II and NEXUS 3.97 models brings to proton 
spectrum with break at 10-15 TeV and power index 7j, 
before break close to that, measured in ATIC-2 experi- 
ment. Nevertheless one can see, that absolute proton flux 
for QGSJET 01 is hardly compatible with any data of di- 
rect experiments, and the break for all these three models 
is more moderate, compared to what can be inferred from 
ATIC-2 data, which though become less definite right 
in the break region. On the other hand SIBYLL 2.1 
allows to reproduce shape of the muon spectrum with 
single power law proton spectrum, which is in reason- 
able agreement with both emulsion chamber and ATIC-2 
data within experimental errors. Further improvement of 
our understanding of the situation, which is of primary 
astrophysical interest, can be achieved via experimental 
study as of muon CR component characteristics with wa- 



ter and ice neutrino telescopes and so of inclusive spectra 
p + A ^ TT , K in fragmentation region. Reduction of 
uncertainties in the latter component with the use of the 
data on primary spectra, hadron and muon components, 
does not look possible, because of 1) poor correlation 
between muon and hadron production mechanisms, 2) 
ambiguity of existing CR experimental data and 3) pos- 
sibility to realize self-consistent description of the data on 
muons and hadrons with the models, having remarkably 
differing inclusive spectra and underlying philosophies. 
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